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We present an effective field theory for the Kondo lattice, which can exhibit, in a certain range 
of parameters, a non Fermi liquid paramagnetic phase at the brink of a zero temperature Anti 
Ferromagnetic (AF) transition. The model is derived in a natural way from the bosonic Kondo- 
Heisenberg model, in which the Kondo resonances are seen as true (but damped) Grassmann fields 
in the field theory sense. One loop Renormalization Group (RG) treatment of this model gives a 
phase diagram for the Kondo lattice as a function of Jk where, for Jk < Jc the system shows AF 
order, for Jk > Ji one has the heavy electron phase and for J c < Jk < Ji the formation of the 
Kondo singlets is incomplete, leading to the breakdown of the Landau Fermi liquid theory. 
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In the last ten years, there has been an increasing 
body of experimental evidence showing striking devi- 
ations from conventional Landau Fermi Liquid theory 
(LFL) at some heavy fermion quantum critical points 
(QCP) [1-10]. The specific heat coefficient is seen to di- 
verge at the QCP, showing a logarithmic increase [2,7,9,8] 
as the temperature is decreased, followed for two com- 
pounds (YbPdi2Si2 and CeCoIns tuned to criticality un- 
der magnetic field), by an upturn at lower temperature 
[9-11]. The resistivity is quasi- linear in temperature for 
most compounds [5,10,6], with perfectly linear depen- 
dence for YbRh 2 (Sio.95Geo.o5)2 and CeCus.gAuo.i [1,2,9]. 
NMR and /i-SR studies show that, generically, the sus- 
ceptibility acquires anomalous exponents [5,10,8]. For 
CeCug-^Aua; [8] neutron scattering measurements reveal 
the presence of B/T and uj/T scaling in the dynamic 
spin susceptibility at the QCP. Recent e-SR [12], ther- 
mal expansion [13] as well as Hall effect [14,15], heat 
transport [16] and Nernst effect [17] studies show further 
deviation from the LFL predictions. Last, for both the 
compounds CeCoIns [18] and YbRli2(Sio.95 Geo. 05)2 [19] 
driven to criticality by applying a magnetic field, there 
are indications that deviations from the LFL theory of 
metals might appear over a whole region -instead of a 
point- of the phase diagram, at the vicinity of the QCP. 

These very striking results have inspired several theo- 
retical descriptions. Some require fully anisotropic (2D) 
spin fluctuations [20,21]. Recently the idea of a local 
mode at criticality [1,8,21] has emerged in relation to 
the compounds CeCu6-xAu x [8] and YbRli2Si2 [9,10]. 
There are also some approaches invoking deconfinement 
and "fractionalization" in gauge theories [22], but none 
of the theoretical approaches so far, can account for a 
non Fermi liquid paramagnetic phase, as well as can fit 
more than one or two experimental observations. 

In this paper we present an effective model for the 



Kondo lattice, which exhibits a non Fermi liquid para- 
magnetic phase at the vicinity to a zero temperature 
AF transition or QCP. The starting point is the Kondo- 
Heisenberg lattice model, where we use a Schwinger bo- 
son representation for the spins of the impurities. We call 
this model Bosonic Kondo-Heisenberg model (BKH). In 
this framework, the Kondo bound states are represented 
by spinless Grassmann fields (x , x)i which have no ex- 
pectation value at the mean-field level. On the other 
hand, the bosonic representation of the impurity spins en- 
ables us to get a decent treatment of antiferromagnetism 
at the mean-field level. We observe that, starting from 
high energy and going to the infra-red sector, the Kondo 
bound states (x, x) acquire dynamics, damping and dis- 
persion, since they are coupled to "itinerant" excitations. 
The key ingredient of our effective theory is thus to treat 
the Grassmann Kondo resonances as true fields in the 
field theory sense. We attribute to them some dynamics 
as well as some damping from the start, and couple them 
in a physical way to fields of the BKH model. A Renor- 
malization Group (RG) treatment applied to this model 
has the remarkable property that, both the formation of 
Kondo singlets and the AF fluctuations - which prevent 
Kondo singlets to be formed- appear at the one loop level 
in the form of a logarithmic singularity. That the forma- 
tion of Kondo singlets has a logarithmic signature at the 
one loop level is a well known fact, in the context of the 
Kondo impurity, where it has been observed a long time 
ago [23,24]. The new feature, here, is the presence of a 
logarithm at the one loop, characteristic of the AF fluc- 
tuations. This is the direct consequence of treating the 
resonances as true fields, with some intrinsic dynamics as 
well as damping. We are then in a position to obtain a 
phase diagram for the Kondo lattice, which can exhibits 
non Fermi liquid properties over a finite region of the 
diagram. 



1 



The BKH lattice Hamiltonian 

H = H c + Hk+Hh , 
where H c = ^ekfl^fka , 

H K = J K J2blb ia ,flf ia , 

iaa' 



(1) 



describes the conduction band (H c ), the Kondo coupling 
between local moments and the conduction electrons at 
site i (Hk), and the super-exchange between neighboring 
spins (Hh)- 

The physics of this model depends on the ratio x = 
Jk/Jh- In the Doniach scenario [25] for the Kondo lat- 
tice, when x« 1 the spins anti-ferromagnetically order, 
and as x is increased the AF state undergoes a transi- 
tion towards a heavy electron paramagnct. It is the goal 
of this paper to take a closer look at the nature of this 
transition. 

When wc formulate the BKH model as a functional 
integral, we can decouple the fields as follows 



Xj Xi 



H K ^H' K = Y J [blxlU + h.c. 



(2) 



I A,; 



Jh ' 



where the bond variable in the first term is a Grassmann 
field which doesn't carry a spin and the bond variable 
in the second term has been chosen following an SP(2N) 
decomposition of the interaction [26]. The mean- field 
theory of this model requires an additional constraint 
term on the Schwinger boson [27,28] representation of 
the impurity spin 



H = H C + H K + H 



J2\(blb ia -2S) , 



where S is the spin of the impurity, taken to be 1/2 in 
the case of our concern. 

At the mean-field level, the x-fermions don't acquire 
any expectation value, while the bosonic bound state A 
can condense, leading to an AF phase. After diagonaliz- 
ing the bosonic part of H in a mean-field approximation, 
one gets the following propagator for the spinons 



G(oj,k) 



{Tb ka bl k _ a ) {Tb ka b_ k _ a ) 
(Tbl k _ a bl) (Tbl k _ a b- k .„) 
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M 2 



to 



k L 



iui — A Afc 
Afc — iu> — A 



(3) 



where w is a bosonic Matsubara frequency, A^ = 
A(sinfc 2; +sa\k y +sinfc z ) is an odd function of k and 



A 2 .. If |A| > A the system is in the param- 



agnetic phase and spinons have a mass ujq = VX 2 — A 2 . 
When |A| = A the system is at a second order phase 
transition towards an AF ground state and spinons be- 
come massless (ui n = 0). In what follows we consider that 
the Kondo lattice undergoes an AF transition, such that 
part of the phase diagram is paramagnetic ( loq ^ 0) and 
part is antiferromagnetic (uj = 0) . uoq is the mass of the 
spinons and plays the role of a parameter in the effective 
theory. 

Through coupling to the spinons ( b krT ) and the itiner- 
ant electrons (f ka ) the x-fcrmion Kondo resonances ac- 
quire some dynamics, damping and dispersion. It is nat- 
ural to think that at some energy scale the dynamics will 
become linear in u> -actually a propagator linear in cj 
is needed to regulate the high energy dependence. We 
propose the following effective action, comprising three 
fields 

S = Sf + Sb + S x + Si n t , (4) 
S f = / dud d k^2fl a (iuj - e k ) f krT , 

S b = f diyd d qJ2V^G-\u,q)V b , 

S x = J dujd d pxl (iu> + iT{w)sgn{w) + gl/J K ) X P , 

Sint = 9 J(duj) 2 d d kd d qd d pS p+ci ^Y,( b l4f^ + h - C ) > 



where V b 



is the spinor basis for the mean- 



field description (3) of the AF fluctuations, T(lu) is the 
damping associated with the x-fermions, which we take 
constant in order to simplify computations, but which 
in all generality is taken to be relevant in the infrared 
sector ( for example T(ui) = \u\ a with a < 1), and g 
is the interaction between the three fields. We have re- 
scaled the x fields such that their inverse propagator has 
the dimension of energy. Now the x-fermions are consid- 
ered as true fields with a finite lifetime and a dispersion 
e x (p) = —g^/Jx- The dispersion becomes a tuning pa- 
rameter of our theory. Generically, three cases are pos- 
sible: (1) The x-"band" is full, or J K > 0. This case 
corresponds to AF Kondo coupling and is considered in 
this paper. The dispersion of the band is neglected as 
long as the band is full. (2) The x-"band" is empty, 
or Jk < 0. This is the ferromagnetic Kondo coupling, 
where the interaction g is marginally irrelevant. This case 
is of no great interest for HF compounds, but can be of 
some relevance for the manganitcs. (3) The x-"b an d" 
is so much renormalized through its interaction with the 
spinons and the conduction electron, that it acquires a 
"Fermi" surface-or a surface in fc-space where the disper- 
sion changes sign. This is a very intriguing possibility, 
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allowing the presence of massless xTermions. This case 
will be addressed in a future work. 

We now turn to the study of the stability of the model 
with respect to the coupling g (diagram c in Fig 1). 
One peculiarity of our model is that at one loop, there 
is no renormalization of the vertex coming from irre- 
ducible diagrams. The logarithmic contribution comes 
from the renormalization of the spinon and %-fermion 
propagators (diagrams a and b in figure 1). We use a 



b) 



c) 6 t 




a > X 



FIG. 1. Diagrams associated with a) the renormalization 
of the x-fermion field, b) the spinon field c) the vertex g and 
d) two loops renormalization of the vertex. Solid, dashed and 
wavy lines are respectively /-fermion, x-f erm i° n an< l boson 
propagators. 

RG scheme, where the parameter g 2 /JK as well as A 
are fixed, while the counter-terms are absorbed into the 
renormalization of the fields x = %x X R an d b = Z^ 2 b R , 
where x R ancl b R are the renormalized fields. The in- 
variance of the propagators is ensured by introducing 
two new parameters, keeping track of the field's lifetime 



9qI ' Jk and G b — Z buJ iu — A. The (3- 



function for the local interaction g takes the form 



dg 



dlogD 



-npg 3 {j K /gl - 1/A) 



(5) 



where D is the bandwidth of the conduction electrons, p 
is the density of states and A is generically positive, such 
that for Jh — in (1) spinons have poles at positive 
energies. In the flow (5) the first term comes from the 
renormalization of the Kondo resonances (diagram a in 
Fig. 1) and is responsible for the usual Kondo behavior of 
the model, leading g to strong coupling. The new feature 
is the presence of a logarithm at one loop, coming from 
the renormalization of the spinon propagator (diagram b 
in Fig. 1) and opposing the Kondo effect in the RG flow. 
Note that this term wouldn't be present if we didn't give 
dynamics to the x-fields in (4). Two different possible 
flows are shown in Fig. 2. For JK—g^/X > 0, the flow goes 
to strong coupling. We identify the strong coupling fixed 
point with the heavy Fermi liquid phase, where Kondo 
singlets form. On the other hand, for Jk — g 2 /X < 0, the 



I*<0 




J*>0 




FIG. 2. RG flows for the coupling constant g. 

flow goes to weak coupling and the model is stable. Note 
that the stability criterion doesn't depend on the sign of 

9- 

Integrating (5) enables to get the energy scales: 



!J 



9 2 o 



1 + 2npJ* log(T/D) 



; r*=Z>exp[-l/(27rpJ*)] , 



where J* = Jk — 5o/A, and T* is the temperature at 
which the Kondo singularity appears, or in other words 
the "lattice" Kondo temperature. One notices, that com- 
pared to the single impurity Tk — Dexp [— 1/(2itpJk)], 
the lattice Kondo temperature T* is reduced by the pres- 
ence of AF fluctuations. 

Examination of the field's lifetimes leads to 



dZ 



dlogD 



-2irpJ K g 2 /g 2 



dZ, 



dlogD 



2irpg 2 /X (6) 



The lifetime of the x-fermion becomes infinitely small in 
the IR sector, which signals a break down of the Landau 
Fermi liquid theory, while at one loop the spinon lifetime 
is stable. 

Three phase diagrams are possible, depending on the 
strength of the AF fluctuations. At large Jk we are, at 
zero temperature, in a heavy Fermi liquid phase, where 
our model flows to strong coupling. The first possibility 
is that J* stays positive in the whole paramagnetic phase, 
reaching the AF QCP while the heavy quasi-particles are 
still well formed. In this scenario, conduction electrons 
are destabilized by a Spin Density Wave (SDW) [29] con- 
necting portions of the Fermi surface. The best candidate 
for this scenario is surely CeNi2Ge2, for which a study 
of the Gruneisen parameter has shown reasonable agree- 
ment with 3D Moriya theory [13]. 

The second possibility is that J* stays positive in the 
whole paramagnetic phase, but vanishes at the AF QCP, 
where the fluctuations are the strongest. This happens 
for J K = gl/X m i n = g%/A = gl/{SJ H ), in a large S 
expansion of the Heisenberg model. 

Last, it might happen that the AF fluctuations close 
to the QCP are so strong that the Kondo effect is bro- 
ken already in the paramagnetic phase, when J* changes 
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Phase diagram for Jjc < So /A 



sign before reaching the AF transition (see Fig. 3). This 
happens for Jk < g 2 /^-- In this case, the Lan- 
dau Fermi liquid theory is broken in a whole region 
of the phase diagram for which our model (4) is sta- 
ble. Recently, measurements on CeCoIns [17,18] and 
YbRh2(Sig.95 Geo. 05)2 [19] indicate the presence of such 
a phase. 

In order, to further justify our identification of the 
strong coupling with the heavy Fermi liquid phase, we 
derived the RG flow at two loops. The diagram d in 
figure 1 is the sole contributor, and we find 

0(g) = ~2irpg 3 (J* + 2irpg 2 A \og{g 2 /J K )) , (7) 

where A = 87rA (u) m — loq)/v s ) 2 in D = 3, with w m and 
wo the upper and lower energies of the spinons and v a 
the spinon dispersion uj q = v s q taken linear over the 
whole spectrum. Since at high energy the starting pa- 
rameters of our model are g 2 /J^ > lj the correction 
to scaling has the sign of the AF Kondo coupling Jk- 
When J* > 0, the two loop flow reinforces the tendency 
to flow to strong coupling, and comforts our identification 
of the strong coupling fixed point with the heavy Fermi 
liquid. When J* < two loop corrections have opposite 
sign, revealing the presence of an intermediate energy 
scale Ti, where the f3- function changes sign. When one 
rcnormalizes from high energy, the system first wants to 
form Kondo singlets, till we reach T±. Further down in 
energy, for J > J\ the formation of the Kondo heavy 
quasi-particles continues, while for J < J± it stops and 
the model (4) is stable, where the three fields correspond- 
ing to damped spinless ^-fermions, conduction electrons 
and spinons interact weakly with each other. The energy 
scale T\ is seen in most of heavy Fermion compounds. It 
can be identified as the scale above which Curie suscepti- 
bilities are observed, or the scale at which the resistivity 
and the specific heat coefficient have a maximum. 

The precise determination of the nature of the Non 
Fermi Liquid (NFL) phase deserves more intensive work. 
The reason why so many compounds- with low disorder- 
show linear resistivity over a wide temperature range may 
lie in the conjecture that a damping T(lu) = |w|log(w) 
stabilizes our model at an intermediate energy scale. 



It has been advanced recently [22] that deconfinement 
of spinons and the emergence of a "fractionalized" state 
might be the clue for understanding the heavy Fermions 
phase diagram. We believe, however, that the key most 
probably lies in the idea of a competition between the 
formation of Kondo singlets and AF fluctuations, which 
are maximal close to a QCP. Our study of model (4) 
shows that de-confinenent of spinons [30-32] is a sepa- 
rate issue from Kondo singlet formation. Even though 
our mean-field treatment of the Heisenberg part assumes 
spinon deconfinement, previous work [33] incites us to 
believe, that in D = 3, gauge theories being deconfining, 
spinons are de-confined in the NFL phase. The question 
is less clear in D = 2. 

Since in our treatment, we have assumed that the \~ 
band is full, we have not addressed the question of the 
variation of the Fermi surface volume. In all generality, 
though, the chemical potential of the x-fermion can be 
renormalized to zero, leading to a reconfiguration of the 
charge carriers close to the QCP. 

In conclusion, we have introduced an effective theory 
for Kondo lattices where the Kondo bound states are 
considered as true fields, with the statistics of spinless 
Fermions. RG analysis of this model signals both the 
formation of Kondo singlets and AF fluctuations as loga- 
rithmic singularities at the one loop level. This leads to a 
phase diagram where direct competition between Kondo 
and AF is apparent, resulting in the possibility of a NFL 
phase separating the AF phase from the heavy Fermi liq- 
uid. We believe the nature of the elementary excitations 
at the QCP in heavy Fermion compounds is essentially 
captured by (4). 
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